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(57) Abstract: Described is a transmission system for transmitting a multicarrier signal from a transmitter (10) to a receiver (20). 

^ The multicarrier signal comprises a plurality of subcairiers. The receiver (20) comprises a channel estimator (28) for estimating 
amplitudes of the subcarriers and for estimating time derivatives of the amplitudes. The receiver (20) further comprises an equahzer 
(24) for canceling intercarrier interference included in the received multicarrier signal in dependence on the estimated amplitudes and 
derivatives (29). The receiver (20) comprises a multiplication by a leakage matrix , and wherein the multiplication is implemented 

^ as a sequence of an N-point IFFT (82), N pointwise multiplications (84) and an N-point FFT (86). 
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Multicarrier transmission system with reduced complexity leakage matrix multiplication 



The invention relates to a transmission system for transmitting a multicarrier 
signal from a transmitter to a receiver. 

The invention further relates to a receiver for receiving a multicarrier signal 
from a transmitter. 

5 

Multicanier modulation methods, such as OFDM and MC-CDMA, have been 
around for some time now. OFDM or Orthogonal Frequency Division Multiplexing is a 
modulation method designed in the 1970's in which multiple user symbols are transmitted in 

10 parallel using different subcairiers. These subcarriers have overlapping (sinc-shaped) spectra, 
nonetheless the signal waveforms are orthogonal. Compared to modulation methods such as 
BPSK, QPSK or MSK, OFDM transmits symbols which have a relatively long time duration, 
but a narrow bandwidth. Mostly, OFDM systems are designed such that each subcarrier is 
small enough in bandwidth to experience frequency-flat fading. This also ensures that the 

1 5 subcarriers remain orthogonal when received over a (moderately) frequency selective but 
time-invariant channel. If the OFDM signal is received over a such channel, each subcarrier 
experiences a different attenuation, but no dispersion. 

The above mentioned properties of OFDM avoid the need for a tapped delay 
line equalizer and have been a prime motivation to use OFDM modulation methods in several 

20 standards, such as Digital Audio Broadcasting (DAB), the Digital Terrestrial Television 
Broadcast (DTTB) which is part of the Digital Video Broadcastmg standard (DVB), and 
more recently the wireless local area network standard HIPERLAN/2. Particularly in the 
DAB and DTTB applications, mobile reception under disadvantageous chaimel conditions 
are foreseen, with both frequency and time dispersion. Mobile reception of television has not 

25 been regarded as a major market up to now. Nonetheless, the DVB system promises to 

become a high-speed delivery mechanism for mobile multimedia and internet services. At the 
IFA '99 Consumer Electronics trade show, a consortium of Nokia, Deutsche Telecom and 
ZDF demonstrated mobile web browsing, email access and television viewing over an 
OFDM DVB link, with a GSM retum channel. With 8k OFDM subcarriers, over the air DVB 
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reception functioned properly for vehicle speeds upto 50 mph. Mobile reception, i.e. 
reception over channels with Doppler spreads and the correspondmg time dispersion remains 
one of the problems associated witii OFDM systems in particular and multicarrier 
transmission systems in general. Whereas its robustness against frequency selectivity is seen 

5 as an advantage of OFDM, the time-varying character of the channel is known to pose limits 
to the system performance. Time variations are known to corrupt the orthogonality of the 
OFDM subcarrier waveforms. In such a case, Intercairier Interference (ICI, also referred to as 
FFT leakage) occurs because signal components from one subcarrier cause interference to 
other, mostly neighboring, subcarriers. 

10 In the document "Equalization of FFT-leakage in mobile DVB-P', Master 

Thesis in Radiocommunication from the Royal Institute of Technology, Stockhohn, by 
Guillaume Geslin, April 1998, a multicarrier transmission system is disclosed. In this known 
transmission system ICI is cancelled (i.e. detected and removed from the received 
multicarrier signal) in the receiver by means of an equalizer. This equalizer derives a vector 

1 5 of estimated symbols from a vector of received symbols. The operation of the equalizer is 
based upon a channel model in vMch the amplitudes of the subcarriers and the time 
derivatives thereof are indicative of the ICI. The receiver comprises a channel estunator 
which generates estimates of these amplitudes and derivatives and supplies these estimates to 
the equalizer. The equalizer then cancels the ICI in dependence on the estimates of the 

20 amplitudes and derivatives. The receiver m the known transmission system is relatively 
complex, i.e. a relatively large number of computations is needed to implement the channel 
estimator and the equalizer. 



25 It is an object of the invention to provide a transmission system according to 

the preamble in which the computational burden is substantially reduced. This object is 
achieved in the transmission system according to the invention, said transmission system 
being arranged for transmitting a multicarrier signal from a transmitter to a receiver, the 
multicarrier signal comprising a plurality of subcarriers, the receiver comprising a chaimel 

30 estimator for estimating amplitudes of the subcarriers and for estimating time derivatives of 
the amplitudes, the receiver fiirther comprising an equalizer for canceling intercarrier 
interference included in the received multicarrier signal in dependence on the estimated 
amplitudes and derivatives, wherein the receiver comprises a multiplication by NxN 
leakage matrix H , and wherein the multiplication is implemented as a sequence of an iV- 
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point IFFT, iSTpointwise multiplications and an iV-point FFT. The invention is based upon the 
recognition that the complexity of the multiplication by the leakage matrix E , which is an 
operation that is used quite conunonly in the context of ICI cancellation in multicarrier 
transmission systems, may be reduced due to the fact that the leakage matrix S is 

5 diagonalized by a Fourier basis, i.e. that E = FH^F" , where F is the iV-point FFT matrix with 
normalized columns and A is a positive diagonal matrix. Hence, a multiplication by the 
Ny^N matrix E may be implemented as a sequence of an N-^omX IFFT, N pointwise 
multiplications and an iV-point FFT, thereby substantially reducing complexity. 

In an embodiment of the transmission system according to the invention the 

1 0 receiver is a decision feedback receiver and wherein the channel estimator comprises a 
smoothing filter for smoothing the estimated amplitudes and/or derivatives. Application of 
such a smoothing filter has the advantage that it exploits the correlation among derivatives. 
That is, since an esthnate of a derivative on a particular subcarrier is inaccurate because of 
noise or other effects, it is usefiil take also into account the values of the derivative at 

1 5 neighboring subcarriers. In practice this typically means that one smoothes the values of 
subcarriers of the various subcarriers. 

In a further embodiment of the transmission system according to the invention 
the FFT is further arranged for demodulating the received multicarrier signal. In this way the 
FFT may be used for demodution and as part of the multiplication of the leakage matrix. 

20 The above object and features of the present invention will be more apparent 

firom the following description of the preferred embodiments with reference to the drawings, 
wherein: 



25 Fig. 1 shows a block diagram of a transmission system according to the 

invention. 

Figs. 2 and 3 and 5 and 6 show block diagrams of embodiments of decision 
feedback receivers according to the invention, 

Fig. 4 shows some graphs illustrating the performance of the decision 
30 feedback receiver as shown in Fig. 3. 

Fig. 7 shows a block diagram of a MC-CDMA transmitter. 

Fig. 8 shows a block diagram of an embodiment of a MC-CDMA decision 
feedback receiver according to the invention. 
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In the Figs., identical parts are provided with the same reference numbers. 

The invention is based upon the development of a simple and reliable channel 
representation. In order to do so, we will consider a multicarrier transmission system, e.g. an 
5 OFDM or MC-CDMA transmission system, with N subcarriers spaced by f, . Each subcarrier 
has a rectangular envelope of a finite length that, including the cyclic extension, exceeds 
(1/ /J.Let^ =[5j,...,5^]^beavcctor of iV^ transmitted symbols, then the transmitted 
continuous time baseband signal may be written as follows: 

10 In the case of a frequency selective time-varying additive white Gaussian 

noise (AWGN) channel, the received continuous time signal may be written as follows: 

wherein the coefficient Hj, (t) represents the time-varymg frequency response at the k -th 
subcarrier, for 1 ^ Jt ^ AT , and wherein n(t) is AGWN within the signal bandwidth. We 
1 5 assume that the channel slowly varies so that only a first order variation may be taken into 
account within a single data block duration. In other words, we assume that every (/) is 

accurately approximated by 

H,{t)^H,{t,)^H,{t,)(t-t,\ (3) 

wherein H],{t) is the first order derivative of //^ (0 and wherein is a reference time within 
20 the received data block. Note that the time varying channel //^ (/) may also take into account 
a residual fi^quency offset, after the coarse frequency synchronization. 

The received baseband signal is sampled with a sampling offset and a rate 
Nf, and a block of its N subsequent samples [y{t^\y{t^ +T\...,yit^ (with 

T = — ) is subject to a fast fourier transform (FFT) of size N. L&ty= [y^ ,...,y^ f be the 
25 vector ofN FFT samples so that 

= T7li>'(^ + nT)exp(-i27ikn/N) . (4) 
After substituting (2) into (4) and using the approximation (3), we obtain 
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I' ~ 

a, =exp(^W^o)(^/(^)+^;(^)(^o -^)). (6) 
d,=txp{}2tf.K)TH]{t,), (7) 

wherein , for 1 ^ ifc ^ , are the samples of AWGN having a certain variance cr^ It is 
5 convenient to rewrite the result (5) in a close matrix form. To this end, we define diagonal 
matrices A = diag {a a„},D = diag {d, } and an iV x matrix 

S = {S,.,};.,.,. B.^^=Y^(.nlN)ixv{ri2n{p-q)nlN). (8) 

With this notation, the expression (5) is equivalent to 

y^As-^BDs-k-n, (9) 

10 wherein n = [w, ,...,Ai;y ]^ is an x 1 vector of AWGN. In the channel model (9), the effect of 
the channel is represented by two sets of parameters a- [a, ,.-,^//]^ and d- [<^p— ,^^;/]^- 
Check that + hence the coefficients aj^,forl<^<iV,are 

equal to the complex amplitudes of the channel frequency response rotated by the sampling 
phase exp(i2;^,/fo) . Similarly, the coefficients , fori ^k^N.aiQ equal to the time- 

1 5 domain derivatives of the complex amplitudes of the chaimel frequency response scaled by 
the sampling period T and rotated by the same sampling phase exp(i2ffJtQ ) . 

Note that an inter-carrier interference occurs when the chaimel response varies 
in time (i.e. d^O). This interference is defined by the vector d as well as the fixed NxN 
matrix H . It can be is easily seen that according to (8) the latter matrix is a Toeplitz 

20 Hermitian matrix and that 

p = q; 

Later in this document, we will call a the (vector of) amplitudes, d the (vector 
of) derivatives and H the leakage matrix. In the above expression, the values on the diagonal 
of S depend on the (arbitrary) choice of the reference time instant to, and can therefor vary 
25 depending on the embodiment of receiver. Typical choices for to are time begiiming, the end 
or middle of a frame window. For to chosen near the middle of the frame, the diagonal terms 
tend to become approximately zero. 

For the implementation of a receiver based on the principles discussed here, 
multiplication by S may be prohibitively complicated particularly for large N (many 
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subcarriers). One can of course only use the tenns near the diagonal and exploit the Toeplitz 
character of H by implementing it as a delay-line filter. But there are more efficient 
implementations of H . We note that the first-order ICI terms result from amplitudes linearly 
increasing with time. That is, one can implement S as the cascade of 
^ • an I-FFT operation (to go back from a frequency domain representation of the 

subcarriers to time domain), 

2- a multiplication of the resulting time-domain signal by a diagonal matrix, i.e., 
a weighting of each component by a scalar. In its basic form the weighting is y 
a linearly increasing function, but m practice one may divert from this, for 
instance to jointiy optimize the reduction of ICI and the avoidance of noise 
enhancements. 

3- an FFT operation to return to the frequency domain representation in which 
subcarrier signals are typically processed. 

This allows the implementation of E with a compl&dty NlogN, instead ofl^ multiplications. 
The inverse of S can be implemented with this structure. For tiie inverse the terms on the 
diagonal of tiie matrix of step 2 will approxunately be of the form or n/(c + n% where n is the 
index at the ou^ut number of the I-FFT. We will exploit tius in several receiver 
embodiments to be given later. In a practical implementation, it is useful to use the same 
hardware circuit for the FFT and IFFT used for S as the hardware used for the main FFT 
operation to retrieve subcarrier signals. 

To process the received signal, the set of channel parameters amid should be 
estimated. The estimation accuracy of tiiese 2N scalar parameters may be enhanced if the 
statistical properties of the channel are used. First of all, we assume that channel variations 
are slow enough so tfiat //^ (r) do not change substantially witiiin tiie duration of a symbol. 
In this case, we may rewrite (6) and (7) as follows: 

d, «exp(/2;5rAW(^o), 1^/^^. (10) 
Let us analyze the relationship between the quantities a, d and physical 
parameters of die propagation channel, namely Uie set of its A: propagation delays {ro,...,r^ }, 
tiie corresponding Doppler shifts {/o,...,/, }, and complex amplitudes {h^ } . Note that 
the statistical properties of die channel frequency response depend on tiie relative delays and 
Doppler shifts whereas tiie group delay and/or Doppler shift result in rotations of , for 
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l^k^K, where the rotations are taken care of by time and carrier synchronizationAracking. 
Hence, we may assume without loss of generality that to = 0 and /o = 0 . Now, the channel 
frequency response Hf and its derivative H] may be written as follows: 

llsO 

H'iif)=i2n%fX^x^2n{f„t-f,lz„)),\^l^N. (11) 

11 BO 

The relationships (10) and (11) may be readily used to deduce the statistical 
properties of the amplitudes a and derivatives d. Whenever the number of propagation paths 
is big enough (ideally i:»7\^), the set of coefficients {///(/)^/f;(OL^^ may be considered 
jointly Gaussian distributed. Moreover, one can show that the sets {h^ it%^<j4 and 
fe (OLiAf are mutually uncorrelated when {AJ,^^^ arc uncorrelated and the Doppler 
spectrum has a symmetric shape. In this case, the vectors a and dmsy be assumed 
statistically independent multivariate Gaussian with zero mean and covariance matrices 

B^a"]^C^,^{dd"}^C, (12) 

where E{.} stands for the mathematical expectation operator and C^,Qare NxN Hermitian 
non*negative definite matrices. 

An important particular case of C„,C^ corresponds to a standard model for 
mobile channels, as described in the book Microwave Mobile Communications by C. Jakes, 
John Wiley & Sons, Inc., 1974. This model (known as Jakes model) assumes mdependent 
contributions of different propagation paths, an exponential delay profile and uniformly 
distributed angles of incidence for different paths. One can show that in this case, 

wherein is the magnitude of the Doppler spread and wherein 7; is the root mean square 
propagation delay spread. The last two parameters depend on the mobile velocity and 
propagation environment respectively. 

Although the outlined channel model is characterized by W parameters, the 
number of independent degrees of fi-eedom is substantially smaller in practice. This property 
comes fi-om the fact that the propagation delay spread is often much smaller than the word 
duration. This property also means that the entries of a are strongly correlated, to the extend 
that the covariance matrix may be accurately approximated by a low-rank matrix. 



f 
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Similarly, the entries of are strongly correlated and the covariance matrix may also be 
accurately approximated by a low-rank matrix. Let us consider the Jakes model and therefore 
(13). Define the eigendecomposition of C : 

C = C/A^/^ . (14) 

5 wherein U is the NxN unitary matrix of eigenvectors of C and wherein A is the NxN 
positive diagonal matrix of its eigenvalues {A, A^ } . Assume that the eigenvalues are 
ordered so that sequence of {A, A^^ } is non-increasing. Under Jakes model, the elements 
of this sequence have an exponentially decaying profile: 

A, ^ exp(- fjr^k), fox\^k<N. (15) 

1 0 Hence, the sequence of eigenvalues may be accurately approximated with a relatively small 
number r of non-zero values: 

{Ai,...,A^}« {a„...,A„0.„0}. (16) 
The aforementioned properties of the channel parameters (i.e. amplitudes and 
derivatives) can be extensively used to derive reduced complexity procedures for channel 

1 5 equalization with ICI removal. Evidently, in situations where the statistical channel may 
divert from the idealized theoretical situation, theses models may still inspire the design of 
practical receiver. The mismatch between the actual channel and the idealized channel model 
may lead to a (small) performance penalty. However, this does not mean that the receiver 
principles disclosed in this invention, can not be used successfully. 

20 Fig. 1 shows a block diagram of a transmission system according to the 

invention. The transmission system comprises a transmitter 10 and a receiver 20, The 
transmission system may comprise further transmitters 10 and receivers 20. The transmitter 
10 transmits a multicarrier signal via a wireless channel to the receiver 20. The multicarrier 
signal may be an OFDM signal or a MC-CDMA signal. The receiver 20 comprises a 

25 demodulator 22 for demodulating the received multicarrier signal, which received 

multicarrier signal comprises vectors of received symbols. The demodulator 22 may be 
implemented by means of a FFT. The demodulated multicarrier signal is supplied by the 
demodulator 22 to an equalizer 24. The equalizer 24 cancels intercarrier interference which 
may be included in the received multicanier signal. The equalizer 24 outputs vectors of 

30 estimated symbols 25 (which have been derived from the vectors of received symbols) to a 
(soft) slicer 26. The slicer 26 produces soft metrics (soft decisions) and/or binary estimates 
(hard decisions) of the (coded) bits to be used in the further signal processing parts of the 
receiver (which are not shown), e.g. a FEC decoder. The output signal of the slicer 26 may 
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also be regarded as comprising estimated symbols 27. The receiver 20 further comprises a 
chamiel estimator 28 for estimating amplitudes 29 of the subcarriers and for estimating time 
derivatives 29 of the amplitudes. The equalizer 24 cancels the mtercarrier interference 
included in the received multicarrier signal in dependence on the estimated amplitudes and 
5 derivatives 29 which are supplied by the channel estimator 28 to the equalizer 24. The 
channel estimator 28 may comprise a reduced complexity filter for deriving vectors of the 
estimated amplitudes and derivatives 29 from the vectors of received symbols 23 and vectors 
of estimated symbols 27. 

We now continue with an example embodiment of a receiver based on the 
1 0 developed channel model. If an OFDM receiver is extended such that it can not only reliably 
estimate amplitudes a (as conventional receivers do), but also (complex valued, e.g. 
including phase information) derivatives d (which is not common for normal OFDM 
receivers), then the user data can be recovered as follows: 

• create the matrix Q^a-^Bd , with a and d denoting the estimates of amplitudes and 
15 derivatives, respectively. Note that the receiver receives the signal y = Qs + n (according 

to (9)). 

• Then estimate s as § = Q'y. 

Such a receiver is called a Imear receiver. The receiver 20 as shown in Fig. 1 may be 
regarded as such a linear receiver when the equalizer 24 implements the matrix multiplication 
20 Qy Here Q* plays the role of an inverse of Q. At least two approaches come to mind. In a 
zero-forcing approach, Q* is the strict algebraic inverse of Q. In an MMSE setting Q' is 
chosen to ensure that § - Q'y = Es|y, i.e. the conditional expectation of s given y. This 
typically minimizes the mean square error \s-s f . In a zero-forcing receiver, the ICI is 

effectively cancelled but noise is amplified. This can lead to undesirable results. The MMSE 
25 receiver optimizes the joint strength of noise and residual ICI. This receiver requires an 
adaptive (typically real-time) inversion of a matrix which depends on the instantaneous 

A 

channel characteristics a and d . 

It is also possible to use a so-called decision feedback receiver. The channel 
model presented earlier in this dociunent reveals that one can refme and improve this 
30 decision feedback receiver in several aspects, among them: 

• A feedback loop in which the estimates of the derivatives are more accurate if one 
exploits more knowledge about the statistical behavior of these derivatives, in particular 
correlation. 
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• A feedback loop in which the error correction code is used wdth in the loop. As a side 
note we mention that for multi-carrier CDMA, the spreading code plays an identical role 
as the error correction code. That is, one can place the decoding (whether is occurs as 
CDMA despreading, as error correction decoding, or any other form) within the loop. 
Fig, 2 shows a block diagram of an embodiment of a decision feedback 
receiver. The decision feedback receiver 20 comprises a demodulator 22 for demodulating 
die received multicarrier signal, which received multicarrier signal comprises vectors of 
received symbols. The demodulator 22 may be implemented by means of a FFT. The 
demodulated multicarrier signal is supplied by the demodulator 22 to a subtracter 32. The 
subtracter 32 subtracts an estimation of tiie ICI included in the received multicarrier signal 
from the demodulated multicarrier signal. The resulting 'ICI-free' signal is supplied to an 
equalizer 24 for normal equalization of tiie signal and to a channel estimator 28. The 
equalizer 24 may also comprise a slicer. The equalizer 24 operates in dependence on 
esthnated amplitudes which are supplied to the equalizer 24 by the chaimel estimator 28. The 
output signal of the equalizer 24, comprising vectors of estimated symbols, is supplied to a 
multiplier 3 1 . Furthermore, the output signal of the equalizer 24 is also supplied to further 
signal processing parts of the receiver (which are not shown). The channel estimator 28 
estimates the amplitudes and time-derivatives of the subcarriers. The estimated amplitudes 29 
are supplied to the equalizer 24 and the estimated derivatives 29 arc supplied to the multiplier 
31. The multiplier 3 1 multiplies the estimated derivatives and the estimated data symbols and 
supplies the resulting signal to a filter 30 which implements the leakage matrix S . The 
filtered signal which is an estimate of the ICI is then supplied to the subtracter 32. 

Based on this general scheme another decision feedback receiver can be 
devised as illustrated in Fig. 3. Here, the FFT demodulator is not shown (but is thought to be 
present). The signal patii comprising Yq, Yi, Y2, the slicer 26, the forward error control 
decoding 42 and the channel estimation (whether blmd or pilot based) are similar to the 
design of a conventional, state-of-tiie-art OFDM receiver. In ttie receiver described here, we 
introduced a subtraction of tiie estimated ICI (HZ5), using Y| = Yo - 2Z5. Here Z5 is an 
estimation of tiie modulated derivatives d s. The signal patii Z|, Z2, Z3, Z4 estimates tiie 
derivatives of tiie amplitudes, with Z4 = d . The rationale behind tfie circuit is tiiat Zi 
recovers the ICI, because an estimate of the modulated subcarriers is subtracted fi-om Yo. 
Only noise, ICI and an estimation error remain. Filter 50 is used to estimate the modulated 
derivatives from tiie ICI. It inverts die leakage matrix H, fliough not necessarily as flie inverse 
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in strict mathematical sense. Preferable this is done while avoiding excessive noise 
enhancement or enhancement of estimation eiror. Modulation of the derivatives is removed 
in the step Z2 Z3. Filter 54 exploits the correlation between the subcarrier derivatives, to 
generate a better estimate Z4. A useful implementation of H involves the cascaded use of an 
5 (Inverse) FFT, a multiplication (with a diagonal matrix) and an FFT. 

Although the circuit is depicted as hardware bmlding blocks, a typical 
unplementation may involve iterative software processing. We experimented with an 
iteration method comprised of the following steps for iteration round i: 

• Input: observation Yo, as well as the (i-l)-th estimate of amplitudes fl(i-l), derivatives 
10 d (i-1), and data §(i-l). Here values between brackets denote the number of the iteration 

round. 

• Calculation of Yi(i), using previous estimates of derivatives d (i-1) and data s(i-l), using 
Y,(i) = Y-S(rf(i-l)s(i.l)) 

• New estimate amplitudes a (i) fcom Yi(i), and (not depicted) possibly exploiting 
1 5 knowledge of amplitudes and derivatives in the previous fiame. 

• New estimate of data §(i) 

• Calculation of Zi(i), Z2(i), Z3(i), Z4(i), Z5(i), 

• (Optional step, not depicted in the Fig.) possible use of knowledge of amplitudes and 
derivatives in the previous frame. This step involves operations which exploit correlation 

20 among subcarrier derivatives. 

• Output: new estimate of amplitudes a (i), derivatives d (i), and data s(i) 

Starting condition is the all-zero vector for a (0), d (0) and §(0). 

The filter 50 attempts to recover an estimate of ds from Zi by filtering Z2 = 
MiZi. One mathematical approach is to use the orthogonality principle for an MMSE 

25 estimate. In this case, an appropriate choice for Mi follows from the requirement E[(Z2 - ds) 
Zi"] = On.* We define e as the vector of decision errors, with e = as- a §.This gives 
Ml = {E[rfs(rfs)"]S" + E[^/se"]} [ E[Eds{dsfE^ -^-h a„^ + ee" + Srfse" + (rfs)"s"e] ]'\ 
wherein an is the variance of the noise. Modeling and (pre-)calculating some of the statistical 
expectation values here can be done, but may not be practical for receiver designers. So next 

30 we will search for simplifying approximations. 

One can simplify the resulting M| as Mi = H"[HS" + G ] wherein G is 
empirically determined as G c\In with a constant ct which may be adapted to specific 
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propagation environments, for instance the average BER, the average SNR or the speed of 
the mobile receiver. 

Z3 approximates d , however it contains error contributions due to AWGN and 
estimation error in ^5 and S. Here we can exploit statistical knowledge that we have 

5 developed about the channel behavior, e.g. on correlation of derivatives. The circuit from Z2 . 
multiplication by 1/x, Z3, M2, Z4 to multiplier Z5 is intended to perform this task. The 
multiplications aim at removing and reinserting the data modulated onto the signal. A 
smoothing operation M2 occurs in between. An MMSE filter to estimate Z4 as closely 
approximating d follows from the orthogonality principle E(Z4 - d)Z3 =0/^, thus M2 = 

10 Ed Z3"[EZ3Z3"]'^ . In practice one may find it acceptable to crudely approximate M2 
=E J "[Ed " + R3] ^ Experiments revealed that R3 = C2 In with a constant C2 is a 
workable solution. 

Fig. 4 shows some graphs illustrating the performance of the decision 
feedback receiver as shown in Fig. 3. The strength of the amplitudes and the derivatives (in 

15 dB) is plotted against the subcarrier number. Graph 60 shows the strength of the actual 

amplitudes, while graph 62 shows the strength of the estimated amplitudes. Graph 64 shows 
the strength of the actual derivatives, while graph 66 shows the strength of the estimated 
derivatives. It can be seen that the amplitudes are estimated very well by the decision 
feedback receiver of Fig. 3, while the estimated derivatives deviate somewhat firom the actual 

20 derivatives. 

Fig. 5 shows another decision feedback receiver based on a MMSE structure. 
It allows iterative computing to minimize the variance of the error between input and output 
of the slicer. The slicer box, as used in Fig. 5 may or may not mclude error correction 
decoding. The slicer outputs the esthnate of the data §, as well as 1/S and as . For OFDM 
25 with QAM modulation, l/§ typically is found from a look-up table. Iterative software 
processing can also used in the receiver of Fig. 5. The iteration method contains the 
following steps for iteration round i: 

♦ Input: observation Yo, as well as the (i-l)-th estimate of amplitudes a (i-1), derivatives 
d (i-1), and data s(i-l) 

30 • Calculation of ¥2(1), using previous estimates of derivatives d (i-1) and data §(i-l), using 
Y2(i) = y-2(i?(i-l)s(i-l)) 
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• Newly estimate amplitudes a (i) from Y2(i), and (not depicted) possibly exploiting 
knowledge of amplitudes and derivatives in the previous frame. 

• Newly estimate data §(i), and the corresponding values of l/§(i) (e.g. table look-up), and 
a(i)§(i) 

5 • Calculation of Z6(i), Z7(i), Z8(i), Z9(i). Here Zg© acts as (correction of) the estimate d . 

• Integration of Z9(i) over various rounds, for instance d (i) = ad (i-1) + (l-a)Z9(i) 

• (not depicted) use knowledge of amplitudes and derivatives in the previous frame. 

• Output: new estimate of amplitudes a (i), derivatives d (i), and data s(i) 

Starting condition is the all-zero vector for fl(0), d (0) and §(0). 
1 0 One can take the filters 72 and 76 non-adaptive and identical to Ma = M| and 

M? = M2 of Fig. 3, respectively. A practical value for the integration constant can be a = 0.9, 

so rf(i) = 0.9rf(i-l) + 0.1Z9(i). 

It appears that several implementations refines are possible as are shown in 

Fig. 6. Here, the inverse of S is implemented as a Finite Impulse Response FIR filter. The 
1 5 second filter M7 is implemented as an IIR smoothing filter. Lastly, an FFT - weighting - FFT 

filter is applied to create an estimate of the ICL 

Many fiirther improvements are foreseen: Use of amplitude and derivatives of 

previous frames to better estimate the amplitude and derivative. This can be done either as 

indicated as 'optional step' in the algorithms, or takmg the initial condition of the iteration as 
20 an extrapolation of results from the previous OFDM frame, with a (0) for the new frame 

equals a (final ) plus T d (0), the latter corrected from the duration of any cyclic prefix or 

guard interval. 

The filters in the receiver, in particular Mi, M2, Me, M? may in a practical 
receiver be fixed or be chosen from a library of precomputed values. For instance, the 
25 receiver control system may upon its knowledge of the propagation environment choose 
between optimized settings for stationary reception (in which most of the ICI cancellation is 
switched of), slow mobile reception (some ICI cancellation), or fast mobile reception 
(aggressive ICI cancellation). 

Furthermore, adaptive filters could be used. These can use reliability 
30 information about estimates. This can be achieved by adaptive matrices or identification of 
erasures in the estimates. 

MC-CDMA is an extension of the basic OFDM principle. In 1993, this form 
of Orthogonal Multi-Carrier CDMA was proposed. Basically, it applies OFDM-type of 
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transmission to a mxilti-user synchronous DS-CDMA signal. So it is vulnerable to Doppler. 
As illustrated in Fig. 7, we use the following vector notation. For OFDM, vector s of length N 
carries a 'frame' of user data, with s = [sq, J|, . . . where the elements Sn are user 
symbols. In MC-CDMA, s = Cx, where C is an A^x iV code matrix and x = [xo, JCi, ... xmf 
represents a frame of user data. We will refer to x as JV user signals, without explicitly 
identifying whether or not all symbols come from the same end-user. The *^th column of C 
represents the 'spreading code* of user data stream k, and will be denoted as 
(cilO]y .. Ck[N'l]f. A commonly used special case, which we will also consider here, is 
C = iV"^^ WH// where WHat is the Walsh-Hadamard matrix of size x iST . In that case, 
C = = C" so CC = /a^ wherein iNisihc NxN unit matrix. In another special case, 
namely that of C = Ij^, the MC-CDMA system reduces to OFDM. For ease of analysis, we 
normalize the modulation as E[XiXj]* = 8,/, or equivalently E[xx"] = /y/. Then E[ss^] = 
EC[xx"]C" = CC" = /. 

Fig. 7 illustrates such a MC-CDMA transmitter. Frames are created by a 
serial-to-parallel (S/P) conversion of an incommg stream of data by a serial-to-parallel 
converter 90, applying the code spreading by a spreader 92, an I-FFT 94 and a parallel-to- 
serial conversion 96 with prefix insertion. We will address the transmission of a single frame, 
and assume that interframc interference is avoided by choosing appropriate guard intervals. 
Hence, the elements of vectors 5 and x are constant with time. The frame duration, excluding 
any guard interval is Ts, where cOgTs = 2n, 

If the receiver architectures proposed in the previous sections are used for 
MC-CDMA, basically, the EEC is replaced by the (inverse) code matrix CThe receiver 
depicted in Fig. 8 is an extension of the receiver of Fig. 5 to accommodate MC-CDMA 
instead of OFDM. The difference with the OFDM receiver resides in the box drawn with the 
dashed line. For OFDM this basically contains a gain control and a slicer. For MC-CDMA, 
we added in the code spreading matrix C (which is used in the spreaders 106 and 1 10 and the 
despreader 102). The receiver is designed such that all matrixes are fixed (or taken from a 
library of pre-computed values) and that no divisions are necessaiy. An exception is the 
weight matrix W (used in the filter 100), which has channel adaptive settings. 

It can be shown that in an MMSE setting, W only has non-zero components on 
the diagonal, with 

o„fl„ 4- constant 
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wherein the constant depends on the noise floor. Details of the circuitry to estimate 
amplitudes are not shown in Fig. 8, but can be based on known principles. 

For MC-CDMA, the slicer bases its symbol decisions on energy received fiom 
all subcarriers, thus the reliability of estimates § is much more accurate in subcarriers that are 
5 in a fade. 

The principles of the receivers described above can also be combined with an 
FFT which handles more samples than the usual size FFT. One example is the use of a 
fractionally spaced FFT, another one is the double sized FFT. Moreover, one can even design 
a system that separates components received via amplitudes from those received over 
10 derivatives. 

Although in the above mainly an OFDM transmission system is described, the 
invention is also and equally well applicable to other multicarrier transmission systems such 
as MC-CDMA transmission systems. Large part of the receivers may be implemented by 
means of digital hardware or by means of software which is executed by a digital signal 

1 5 processor or by a general purpose microprocessor. 

The scope of the invention is not limited to the embodiments explicitly 
disclosed. The invention is embodied in each new characteristic and each combination of 
characteristics. Any reference signs do not limit the scope of the claims. The word 
"comprising" does not exclude the presence of other elements or steps than those listed in a 

20 claim. Use of the word "a" or "an" preceding an element does not exclude the presence of a 
plurality of such elements. 
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CLAIMS: 



1 A transmission system for transmitting a multicarrier signal from a transmitter 

(10) to a receiver (20), the multicarrier signal comprising a plurality of subcarriers, the 
receiver (20) comprising a channel estimator (28) for estimating amplitudes of the subcarriers 
and for estimating time derivatives of the amplitudes, the receiver (20) further comprising an 
5 equalizer (24) for canceling intercarrier interference included in the received multicarrier 
signal in dependence on the estimated amplitudes and derivatives (29), wherein the receiver 
(20) comprises a multiplication by NxN leakage matrix S , and wherein the multiplication 
is implemented as a sequence of ani^-point IFFT (82), //pointwise multiplications (84) and 
anA^-pointFFT(86). 

2. The transmission system according to claim 1 , wherein the receiver (20) is a 

decision feedback receiver and wherein the channel estimator (28) comprises a smoothing 
filter (76) for smoothing the estimated amplitudes and/or derivatives. 

15 3. The transmission system according to claim 1 or 2, wherein the FFT (86) is 

fiuther arranged for demodulating the received multicarrier signal. 

4. The transmission system according to any one of the preceding claims, 
wherein the multicarrier signal is an OFDM signal. 

20 

5. The transmission system according to any one of the claims 1 to 3, wherein the 
multicarrier signal is a MC-CDMA signal. 

6. A receiver (20) for receiving a multicarrier signal from a transmitter (10), the 
25 multicarrier signal comprising a plurality of subcanriers, the receiver (20) comprising a 

channel estimator (28) for estimating amplitudes of the subcarriers and for estimating time 
derivatives of the amplitudes, the receiver (20) fiirther comprising an equalizer (24) for 
canceling intercarrier interference mcluded in the received multicarrier signal in dependence 
on the esthnated amplitudes and derivatives (29), wherein the receiver (20) comprises a 
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multiplication by NxN leakage matrix S , and wherein the multiplication is implemented 
as a sequence of an iV^pomt IFFT (82), N pointwise multiplications (84) and an JV-point FFT 
(86). 

5 7^ The receiver (20) according to claim 6, wherein the receiver (20) is a decision 

feedback receiver and wherein the channel estimator (28) comprises a smoothing filter (76) 
for smoothing the estimated amplitudes and/or derivatives. 

8. The receiver (20) according to claim 6 or 7, wherein the FFT (86) is further 
1 0 arranged for demodulating the received multicarrier signal. 

9, The receiver (20) according to any one of the claims 6 to 8, wherein the 
multicarrier signal is an OFDM signal. 

15 10. The receiver (20) according to any one of the claims 6 to 8, wherein the 

multicarrier signal is a MC-CDMA signal. 
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